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Abstract. Nuclear fusion cross-sections considerably higher than corresponding 
theoretical predictions are observed in low-energy experiments with metal matrix 
targets and accelerated deuteron beams. The cross-section increment is significantly 
higher for liquid than for solid targets. We propose that the same two-body correlation 
entropy used in evaluating the metal melting entropy explains the large liquid-solid 
difference of the effective screening potential that parameterizes the cross-section 
increment. This approach is applied to the specific case of the 6 Li(d, a) 4 He reaction, 
whose measured screening potential liquid-solid difference is (235 ± 63) eV. Cross 
sections in the two metals with the highest two-body correlation entropy (In and Hg) 
has not been measured yet: increments of the cross sections in liquid relative to the 
ones in solid metals are estimated with the same procedure. 
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The observed pronounced enhancement of the rate of fusion of deuterons implemented in 
metal matrices with accelerated deuterons is well established since the early and almost 
contemporaneous works by the Kasagi group in Sendai urn G], by the Rolfs group in 
Gran Sasso and Bochum [HE] EE0 M, E] and by the Berlin group [10, LI • 12] . Analogous 
enhancements has been seen in fusion between accelerated deuterons on light nuclei in 
the target [13, ITT] . 

The understanding of these enhancements of nuclear fusion rates is necessary for the 
correct evaluation of rates of stellar fusion processes, of plasma fusion in laboratories and 
of low-energy cross sections for nuclei in vacuum. The exact quantitative link between 
cross sections in different environments (stellar, laboratory, other plasmas, vacuum) is 
still under experimental and theoretical scrutiny ca. 

A standard approach parameterizes these enhancements of the rates adding 
a constant negative effective potential to the total energy of the fusing particles. 
Such effective potentials are much higher than the screening potentials in adiabatic 
approximation, which are several tens of eV [16], while Kasagi et al. [13], [T4j report 
rates of fusion of deuterons on lithium targets in terms of effective potentials that are 
several hundreds of eV larger than adiabatic values. The same group has been recently 
measuring rates in a liquid indium matrix, finding a similar behavior m 

These potentials can be interpreted as an additional screening, in agreement with 
the fact that these solid and liquid metals are dense cold plasmas, where the Debye- 
Huckel model should be modified [18] T9] . 

Several models appear to reproduce only partially such experimental results. We 
have suggested explanations m of the effective potential for fusion in solid metal 
matrices based on: 

1) modifications of the Debye-Hiickel approach by spatial fluctuations of the Debye- 
Hiickel radius; 

2) a critical review of the stopping power; 

3) quantum uncertainty effect which broadens the energy-momentum relation and 
increases the effective kinetic energy of the colliding particles following the Galitzskii- 
Yakimets quantum description [21, 221; GSJ. 

Furthermore, we have more recently shown the importance of the correlation 
entropy term in metal melting m ■ 

In the framework where the environment effect is parameterized by adding constant 
potential to the particles in the colliding beam, the negative potential is found [13, [14] 
about 200 — 300 eV larger in the liquid than in the solid metal. In this paper we show 
that this difference can be explained by means of the two-body correlation entropy term 
of the ions or atoms in liquid phase for the specific case of lithium. We suggest that 
this two-body entropy correlation is important and in same case dominant in fusion 
processes in environments different from vacuum. Being the correlation entropy term of 
In and Hg the largest among the metals [25] [26], we evaluate the liquid-solid difference 


The role of correlation entropy in nuclear fusion 


3 


for these two metals, quantity useful for future experiments. 

In Section 2, we give an overview of the experimental situation and outline the 
standard approaches used to explain the experimental results. In Section 3, we report 
our interpretation of the difference between the screening potentials in solid and liquid 
lithium and, furthermore, we evaluate the predicted difference for indium and mercury 
and conclude. 

2. The experimental situation and interpretation of the spread 

In the very-low energy experiments of nuclear fusion in lithium targets, when the center 
of mass energy becomes smaller than 30 — 40 keV, the cross section is sensibly higher 
than theoretical vacuum predictions. The screening potential that parameterizes this 
enhancement depends on the state of the metal target: 

1) in gaseous targets, one needs screening potentials larger than the electron 
screening predicted in the adiabatic limit (20 — 30 eV). In fact the total screening 
potentials measured by several groups for different gases are in the range 200 — 
300 eV [131127]; 

2) in the few experiments with liquid target, the anomalous increase of the screening 
potentials relative to atomic adiabatic screening plus conduction electron energy in the 
Thomas-Fermi approximation is even larger. The total screening potential for lithium 
is rs (520 ± 190) eV against a prediction of about 71 eV (12, H]; 

3) in the many experiments with solid targets, effective screening potentials range 
from 300 to 1300 eV, again values much larger than the theoretical predictions from 
atomic and conduction electrons [28. 29j 130, 31J [32j, 33]. The specific value reported for 
lithium is ~ (310 ± 70) eV [14]. 

4) Kasagi has very recently found [17] anomalous energy spectra of the outgoing p 
and t in the reaction d (d, p) t in liquid Indium. 

In summary the atomic electron screening is not enough to explain the 
enhancement, if its value is deduced following the standard approach: adiabatic or 
sudden approximation and standard Debye-Hiickel model [16] . It appears that target 
deuterons are not at rest, but they have a finite momentum, and that reactions are not 
simple two-body reactions. Kasagi et al. m suggest that results could be explained by 
assuming cooperative colliding processes; following this hint, we plan to study possible 
connections with the quantum momentum-energy uncertainty [2Tj. 

Attempts to explain the experimental results in Sendai PEI E], in Bochum and 
Gran Sasso [2 [5, 6] 0 [8, 9], and in Berlin [10] [12], FL] are based on: 

ionic Debye screening, assuming liquid lithium is a low-temperature high-density 
plasma (Toriyabe et al. [31]); a simplified model of the classical quasi-free electrons with 
an electron screening distance of the order of the Debye length, but Debye screening 
is a cooperative effect of a large number of electrons inside the Debye sphere and this 
model needs a number of electrons equal or smaller than one (Rolfs et al. [35] 22, 9]); 
use of dynamical screening as from [38] (Dappen et al. [32]); deuteron dynamics with 
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migration of electrons from the host metallic atoms to the deuterium ones (Huke at al. 

mm- 

Other authors [2011221139] point out that, as a result of frequent collisions of particles 
in environments like metal matrices or dense plasmas, the complete correspondence 
between total and kinetic energy disappears and the distribution function depends both 
on total and kinetic energy with visible effect on its tail. Momentum distribution is 
fundamental for the fusion reaction rates and deviations from Maxwellian distribution 
play a central role. Therefore two effects have been examined: the increase of 
the deuteron momentum distribution tail due to the quantum energy-momentum 
uncertainty effect and the spatial fluctuations of the Debye-Hiickel radius, the two effects 
lead to larger reaction rates at energies below a few keV. 

In liquid indium bombarded with deuteron molecular beam, where anomalous 
spectra have been detected, Kasagi et al. m introduce a cooperative colliding process. 

In the following we shall focus on the difference of screening potential for the 
6 Li(d, a)' 1 He reaction in liquid lithium compared to the same reaction in solid (or 
atomic/molecular) lithium. Note that most of the large uncertainty on the two potentials 
is a consequence of the large uncertainty on the bare cross section relative to which the 
screening enhancement is defined. For instance, the values for liquid and solid lithium 
targets from Ref. ||I4] are 520 ± 190 eV and 310 ± 70 eV: the difference would naively 
be 210 ± 202 eV. But the uncertainty on the bare cross section cancels in the difference 
and the percentage error is much smaller, about 35%: for instance Fang at al. [13] quote 
235 ± 65 eV for the difference between the liquid and atomic/molecular phase. As a 
reference typical value we use 200 ± 70 eV. 

3. Our approach, results and conclusions 

Liquid metals are characterized by a finite microscopic correlation length L c : spatial 
correlations or long-range order vanish for distances much greater than L c . The non¬ 
trivial interaction of each particle with the environment is very different from the one 
of a quasi-particle in gas: a better description is a particle that moves within a cage 
formed by its neighbors; dimensions of this cage are of order L c (see Wallace [26] 40] for 
a detailed discussion). 

We propose that the cooperative contribution of the Wage atoms in this cage can 
explain the additive ~ 200 eV effective screening potential in liquid lithium relative to 
solid lithium needed to describe the corresponding increase of the fusion cross section. 
The correlation energy of the cage increases the effective energy of the colliding particles: 
this larger energy can be equivalently interpreted as a reduction of the stopping power 
(correlation reduces the transfer of energy to the environment) or as a direct contribution 
of the cage to the relative kinetic energy of the colliding particles. 

In liquid metals the entropy per atom can be written as 

S = S t + S 2 + S x + S Q + S‘ , 
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where Si is the one-body (kinetic energy) term, S 2 and S x are the two- and many-body 
correlation entropies, S® is the quantum contribution and S e is due to the electron 
cloud. Only the one-body Si and two-body correlation S 2 entropies are large and need 
to be considered for our scope. In fact, Wallace [25] 00] 01] has calculated S 2 for many 
liquid metals as the difference between the experimental value of the total entropy and 
the one-body term Si, disregarding corrections from the smaller terms. In turns, the 
difference between the one-body entropy in the two phases is easily determined [25] and 
it cannot influence the two-body reaction rate: only the two-body correlation entropy 
S 2 needs to be considered to explain the different fusion rates in liquid relative to solid 
metals. 

Our model is the following. When a deuteron enters the metal (solid or liquid) 
suffers the stopping power process. In the solid, at the end of its path, the kinetic 
energy of the deuteron is reduced by the stopping power process and a contribution is 
gained from the atomic electron screening in the lithium environment. The gain can 
be calculated by using the modified Debye-Hiickel approach. Furthermore the deuteron 
collision frequency is increased because of the energy-momentum uncertainty (Galitski- 
Yakimets quantum effect) [20]. Solid lithium has the ions localized in the crystal; the 
entropy of the crystal is practically the sum of one-body contribution Si (kinetic energy) 
plus very minor contributions from many-body correlations, from atomic electron cloud 
and from quantum effects. Crystal neither receives or gives energy to the incoming 
deuteron or to the environment except for the contribution given to deuteron from 
atomic electron cloud (contribution calculated by the adiabatic approximation), the 
ionic screening being negligible (in fact since the mass of ions is greater than electron 
mass, positive ions have slow mobility and cannot respond quickly to change [42; 143]). 
When, in the other case, deuteron enters the liquid metal, its interaction with the 
environment is different for two reasons: it travels within a cage of ions/atoms and 
receives an amount of energy (that we can consider equivalent to a screening energy) 
because the ions in liquid metal, although non localized, are correlated. In the liquid 
metal the negative contribution to entropy of the two-body correlation entropy restores 
in part a sort of order and the cage can give an amount of energy to the deuteron at 
the end of the stopping power process or during the travel. The internal energy of the 
cage increases deuteron kinetic energy. In summary: deuteron enters the metal sample, 
travels into the matter submitted to the stopping power process, gains energy from 
the atomic electron (adiabatic approximation and modified Debye-Hiickel approach), 
increases collision frequency by quantum energy-moment uni uncertainty effect and 
furthermore when travels in the liquid metal collects the amount of energy provided by 
its cage during or at the end of its itinerary. This energy cannot be given to deuteron 
by the solid metal, is responsible of the different behavior of Li+d fusion in solid and 
in liquid metal and represents the difference of the internal energy between liquid and 
solid metal of the cage bunch as quantitatively shown below. 

From the energy balance, by calling A E the energy gained by the deuteron and W 
the energy interaction contributions (besides the standard atomic electron screening) 


The role of correlation entropy in nuclear fusion 


6 


due to interactions of d, Li and Be (the d-Li nuclear compound system just before 
fusion) with the environment, we can set 

(W'p - WT 1 ) + (wff - + (W'S - Wb°«‘) -ae = o. 

While the two differences of the second and third parenthesis are negligible together 
VF| o1 , only is non negligible, because of the d-cage interaction that assures an 

amount of energy from the negative two-body correlation entropy term. 

Since the free energy is the same in the two phases (solid and liquid), the entropy 
change must be matched by a change of internal energy: 

A U = f/ liq - U so1 = Wage = AS = Wage k B T m S 2 , (1) 

where we have ascribed the whole entropy difference to the correlation of the Wage 
particles in the cage and the two-body correlation entropy S 2 is in units of k B . When 
it is known, it is better to use the actual temperature of the experiment T instead of 
melting temperature T m < T. It would be interesting to test the linear dependence of 
AU with the temperature. Note that the energy lost locally by the cage is so small that 
it does not modify the macroscopic properties of the liquid. 

Therefore A E = —AU is the additional energy available for the colliding particles 
in liquid 


E — W + Wage — W + Wage W T m S 2 , 


where E k is the kinetic energy of the deuteron inside of the metal reduced by the stopping 
power, whose value can be calculated following the approach explained in detail by 
Coraddu et al [20J. Since the two-body correlation melting entropy is always negative 
(correlations give order to the system) in a liquid metal, W = E + AE = E — Wage > E. 

In lithium Wallace [25] finds that S 2 (Li) = —2.48; given an experimental 
temperature T = 520 K and a difference of effective screening potentials of about 
200 eV, one finds 


N = 

1 y cage 


A E 


1finn A E 2.48 520 K 

1800 x ———— x —— x 


-S 2 k B T m 200 eV -S 2 T 

In addition we know a theoretical expression for S 2 m 


( 2 ) 


W = ~p f 9l(r) In g(r)dr » 

CCLQ 6 


(3) 


where g(r) is the two-body correlation function, Wage is the volume occupied by the 
Wage atoms of the cage and p is the density of the metal in particles per unit volume, 
that for liquid lithium is about 4.4 • 10 22 cm -3 . From Eq. ([2D the volume of the cage can 
be estimated 


i/ = 

v cage 1y rare 


2AE 


P 


y cage 


2.0 • 10" 19 cm 3 x 


pksTm 

A E 


520 K 


200 eV 


x 


T 


(4) 


Alternatively, if we have a way to calculate/independently measure the volume of 
the cage or the correlation length, we could make predictions and/or tests of the model. 
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For instance, Ichimaru |44j considers liquid metals as strongly coupled plasmas with 
mean interaction energy per particle larger than the kinetic energy: the corresponding 
plasma parameter is much greater than one. Using Monte Carlo methods Ichimaru 
determines terms of a power series expansion of the two-body correlation. 

From the theoretical point of view it can be interesting to express the above results 
in terms of x = L cage /a, the dimensionless correlation length in units of the atomic 
diameter a [45] . Expressing the cage of volume as Wage = Z/ 3 age = x 3 a 3 and using 
Eqs. ljlil . (ED, and (J4J) , the difference of potentials 



(5) 


and the number of atoms in the cage 



( 6 ) 


become both functions of the dimensionless parameter x. 

From this last Eq. (ED using pa 3 ~ 1, N CSige ~ 1800 and — S 2 = 2.48, one finds 
x ~ 20. Under the hypothesis that the entire screening is due to the correlation entropy, 
the cage correlation length is about twenty times the diameter of the atom. 

If we want to compare and make estimates for other metals, it is useful to write 
Eqs. (ED and (ED as ratios of the expressions for metal a and for metal b: 


(AU) a _ (iV cage ) a (T) a (S 2 ) a _ (pa 3 ) a (T) a (z 3 ) 


(7) 


V ^ 7 Cl V Cdtifc; 7 Cl \ 7 Cl \ ^ Z / n 7 Cl \ 7 Cl \' Ay 7 Cl 

(A u) h = (A cag e) b (T) b (S 2 ) b = (pa 3 ) b (T) b (x 3 ) b 


Among the liquid metals, liquid indium and liquid mercury have the highest value 
of S 2 : S 2 (In) = —3.39 and S 2 (Hg ) = —3.97 with increments of +37% and +60%, 
respectively, compared to S 2 (Li ) = —2.48. The corresponding values for per 3 are: 
pa 3 (In) = 1.37 and pa 3 (Hg ) = 1.12 with increments of +12% and —10%, respectively, 
compared to pa 3 (Li) = 1.25. 

We do not have the experimental values of A U for indium and mercury. At the same 
temperature of the lithium experiments T = 520 Ii, assuming that the dimensionless 
correlation length x be approximately the same, Eq. (ED yields values of A U about 
10% larger (smaller) for indium (mercury) relative to the lithium value. If it is assumed 
approximatively equal number of particle in the cage Eq. (ED predicts values of A U much 
larger, about 40% (60%) larger for indium (mercury) relative to the lithium value. In 
addition Eq. ([7]) predicts that A U is proportional to the experimental temperature. 

In conclusion we have shown that it is possible to explain the difference between the 
screening energy to be used in solid and in liquid lithium A U by means of the two-body 
correlation entropy term S 2 and we have estimated the possible values of A U for indium 
and mercury, the two metals with the largest two-body correlation entropy. In addition 
we predict that A U scales linearly with the temperature. These evaluations could be 
useful for future experiments as the ones performed recently by Kasagi group. 
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